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Quantum Dynamics of Spin Wave Propagation Through Domain Walls 
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Through numerical solution of the time-dependent Schrodinger equation, we demonstrate 
that magnetic chains with uniaxial anisotropy support stable structures, separating ferromag- 
netic domains of opposite magnetization. These structures, domain walls in a quantum system, 
are shown to remain stable if they interact with a spin wave. We find that a domain wall 
transmits the longitudinal component of the spin excitations only. Our results suggests that 
continuous, classical spin models described by LLG equation cannot be used to describe spin 
wave-domain wall interaction in microscopic magnetic systems. 
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1. Introduction 

Wave propagation in onc-dimcnsional magnets 
through a magnetic domain wall (DW) is an interesting 
topic in quantum many-body physics. A DW separates 
two regions with opposite magnetization. The DW in 
mesoscopic wire can be considered to be self-assembled 
stable nanostructures which is treated as a kind of 
soliton in a continuous medium. Such structure can 
be created or annihilated by some external action.^ 
The manipulation of DW in stripes has already been 
proposed as a way of storing information or even 
performing logic functions, and to offer new types of 
electronics devices^ in which the DW motion carries the 
information along a magnetic wire of submicrometer 
width, with DW velocities up to thousand kilometers per 
second.^ Recently, a direct observation of the pendulum 
dynamics of a DW has been reported.'^ The DW as a 
topological particle has a very small but finite mass of 
6.6 X 10-23kg.3.4 

The structure of DWs and also wave propagation in 
one-dimensional classical spin systems has been studied 
in Refs.^^^ Recently, the interaction between DWs and 
spin waves has attracted a lot of interest. Hertel et al 
showed that the DW induced phase shift of spin waves in 
the Landau-Lifshitz-Gilbert (LLG) model of thin, narrow 
strips, is a characteristic property of such systems.* The 
value of the phase shift of spin waves passing through a 
DW was found to be proportional to the angle by which 
the magnetization of DW rotates in the film plane.* This 
effect might be used as a concept for a new generation of 
nonvolatile memory storage and logical devices.* 

On the other hand, recent progress in synthesizing 
materials containing ferromagnetic chains^^^^ opens new 
possibilities to study the interaction between a spin wave 
and a DW in a microscopic spin chain. Furthermore, 
quantum spin models provide a playground to investigate 
how quantum information can be transferred in quan- 
tum spin networks. "'^^"^^ But, in contrast to the nanoscalc 
phenomena mentioned earlier, on the atomic level, the 
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Fig. 1. Left pictures (a,b,c,d): Spin configurations at time t/r = 
0; Right pictures (e,f,g): Dynamically stable spin configurations 
for the Heisenberg-Ising model (A = 2), at time t/r = 100. (e): 
Ferromagnetic state of the spin chain; (f): State containing a 
DW at the center of the chain; (g): State containing two DWs. 
Right picture (h): Spin configuration at time t/r = 100 for the 
Heisenberg model (A = 1), illustrating the instability of the ini- 
tial domain wall state (d). 



spin dynamics is purely quantum mechanical and in 
such strongly quantum fluctuating systems it must be 
described by the time-dependent Schrodinger equation 
(TDSE). Then, it is of considerable interest to compare 
the properties of spin wave propagation through a mag- 
netic domain boundary in a single spin chain with the 
dynamics obtained in mesoscopic system, in which the 
magnetization is regarded as a classical, continuous vari- 
able. In the nanoscale regime, the DW is defined as the 
boundary of regions with opposite magnetization. On the 
atomic level, a DW may be defined as a structure that is 
dynamically stable under quantum mechanical motion, 
the existence of which has to be confirmed. 

2. Dynamically Stable Domain Walls 

In this paper, we study the stability of DWs and the 
effects of DWs on the spin wave propagation in a chain 
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of N sites on which we place S — 1/2 spins. We solve 
the TDSE to compute the time-evolution of the magne- 
tization at each lattice site. The Hamiltonian of the spin 
chain is given by^^ 
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where the exchange integrals J > and AJ determine 
the strength of the interaction between the x, y and z 
components of spin 1/2 operators S„ = {S^, , S^). We 
solve the TDSE by the Chebyshev polynomial algorithm 
which is known to yield extremely accurate solutions of 
the TDSE, independent of the time step used.^''"^" We 
display the results at time intervals of r = 7r/5J. We 
present results for systems containing = 26 spins only. 
We checked that simulations for = 20 spins (data not 
shown) yield qualitatively similar results. In our numer- 
ical work, we use units such that h — I and J = 1. 

First, we study the stability of DWs. The left panel 
of Fig. 1 shows the spin configurations that we take as 
the initial state {t/r = 0) in the simulation. All the 
results shown in this Letter have been obtained using 
open boundary conditions. We let the system evolve in 
time according to the Hamiltonian Eq.(l) for a long time 
(t/r = 100) and find that the motion generates a dynam- 
ically stable state with DW(s) (see Fig. l(f,g)). The DW 
is defined as the boundary between regions of different 
magnetization but it is not trivial that these boundaries 
exists in the presence of strong quantum fluctuations.^ 
Figure l(e,f,g) shows the dynamically stable spin con- 
figurations obtained by starting from the corresponding 
configuration (a,b,c). 

Whether or not quantum fluctuations destroy the 
DW(s) depends on the value of the anisotropy A. For 
the model Eq.(l), it is well know^^ that quantum fluctu- 
ations destroy the long range order of the ground state 
if — 1 < A < 1 (XY-like) or A = —1 (Heisenberg an- 
tiferromagnet). For A > 1 (Ising-like) the ground state 
exhibits long range order. This property is reflected in 
the stability of configurations that contain one or more 
DWs, except for A = 1. Our numerical simulations show 
that configurations with a DW are dynamically stable if 
A > 1. For comparison, in Fig. 1 we include the case 
A = 1, where the initial DW structure (Fig. 1(d)) is 
destroyed (Fig. 1(h)). Not surprisingly, the destructive 
effect of quantum fiuctuations can be suppressed by in- 
creasing A. Having studied systems with different values 
of A, we found that A = 2 is representative for the quan- 
titative behavior of the anisotropic systems. Therefore, 
in this paper, we present results for A = 2 only. We also 
checked the effect of the boundary condition. We found 
almost the same stable DW structures in the case of pe- 
riodic boundary conditions (results not shown). 

3. Spin Wave Propagation 

We use the configuration at t/r = 100 as the initial 
configuration to study the spin wave dynamics. We gen- 
erate a spin wave excitation by rotating the left most 
spin Si in Fig. l(e,f,g). For reference, we also consider 
the dynamics of the ferromagnet (see Fig. l(a,e)). This 




Fig. 2. (color online) Time evolution of the magnetization 
{S^ (t)) of the Ising-like spin chain with A = 2. The initial config- 
uration (t/r = 0) of each panel (a,b,c) is shown in Fig. l(a,b,c), 
respectively. At the time t/r = 100, the first spin (n = 1) is 
flipped, generating a longitudinal spin wave. 



(a) No DW 



(b) One DW 



(c) Two DWs 




Fig. 3. (color online) Time evolution of the transverse component 
{S^ {t)) of the magnetization, for the same cases as those shown 
in Fig. 2 except that the first spin is rotated by 7r/2 about the 
i/-axis (instead of flip) at the time t/r = 100. 



case without DW can be analyzed analytically, so that 
it also gives check of precision of numerical calculation. 
Actually, we found very small difference between the an- 
alytical results and numerical ones. 

In Fig. 2(a) we show the time evolution of{{S^ (t))} for 
n = 1, . . . N after flipping Si, in the case of the uniform 
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Fig. 4. (color online) Time evolution of the magnetization 
{Sff {t)). (a) spin wave generated by flipping the first spin Si; 
(b) spin wave generated by rotating the first spin Si by 7r/2 
about the y-axis. Solid (black) line: No DW, corresponding to 
the spin configuration Fig. 1(e); Dashed (red) line: one DW, cor- 
responding to the spin configuration Fig. 1(f). Because {Sff (t)) 
is negative in this case, we plot the absolute value to facilitate 
the comparison; Dotted (blue) line: two DWs, corresponding to 
the spin configuration Fig. 1(g). Comparison of (a) and (b) shows 
that the times at which the {Sf^ (t)) reaches one of the minima 
does not depend on method by which the spin wave is generated. 



chain. The thiie evolution of {{S^ (t))} for n = 1,...N 
in the chain with one DW at n — 13, 14 is depicted 
in Fig. 2(b), and Fig. 2(c) shows the results for the 
chain with one DW at n = 10, 11 and another DW at 
n = 17, 18. Hence, we demonstrate that even in the pres- 
ence of a spin wave, the DW structure remain stable. In 
the model Eq.(l), the magnetization in the z-direction 
is a conserved quantity. Hence, by flipping one or more 
spins we change the total magnetization of the initial 
state. The expectation value of the transverse spin com- 
ponents is identically zero {{S^ (t)) — {S^ {t)) = for 
n = l,...N), for all t > 0. 

From Fig. 2, we can deduce how the spin wave is scat- 
tered by the DW(s). The triangular pattern in Fig. 2(a) 
merely results from the reflection of the spin flip excita- 
tion by the other edge of the chain. The triangular pat- 
tern is also present in Fig. 2(b), but the presence of the 
DW causes (S^ (t)) to change sign if n > N/2. Fig. 2(b) 
also demonstrates that the DW itself is extremely robust, 
even in systems with one spin flipped. A similar behav- 
ior is observed for the case of two DWs (see Fig. 2(c)), 
indicating that the change of sign at the DW is generic. 

The slope of the line in Fig. 2 from the point {n — 
1,</t = 100) that connects spin 1 and spin N is directly 
related to the velocity of the excitation. We can estimate 
the time of the excitation to propagate from site n to site 
m by analyzing the infinitely long chain. Starting from 
an initial state in which we flip the spin at site n, the 
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Fig. 5. (color online) Time evolution of (S^ (t)) of the same sys- 
tem as in Fig. 3, plotted on two different scales. Solid (black) 
line: No DW, corresponding to the spin configuration Fig. 1(e); 
Dashed (red) line: One DW, corresponding to the spin configu- 
ration Fig. 1(f). Dotted (blue) line: Two DWs, corresponding to 
the spin configuration Fig. 1(g). 



magnetization at site m is given by 
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where |m) denotes the ferromagnet state with a flipped 
spin at site m and Jm{x) is the Bessel function of the first 
kind of order m. Although Eq. (3) is valid for the infinite 
chain only, we may expect that it provides a qualitatively 
correct description of the wave propagation in the finite 
system. Our numerical calculations (results not shown) 
demonstrate that for N > 16, the time for the excitation 
to travel from n = 1 to m = 26 agrees within 2% with 
the first minimum of Eq. (3). 

Although it is clear that the longitudinal motion of 
the spin that results from the spin flip can easily prop- 
agate through the DW structures, quantum fluctuations 
reduce the amplitude of the excitation and for t/r > 250 
it becomes difficult to follow the excitation in Fig. 2(b,c). 
As mentioned earlier, we could increase A to reduce the 
quantum fluctuations but this does not change the qual- 
itative features that we are interested in. 

Next, we study the propagation of the transverse com- 
ponents, that is the x or y components of the expectation 
values of the spins. At t/r = 100, we excite the system 
by rotating the first spin in Fig. l(e,f,g) by 7r/2 about 
the y-axis. After this rotation, the magnetization of spin 
Si is parallel to the x-axis. Starting from this configura- 
tion, the time evolution will cause the first spin to rotate 
about the z-axis (due to the presence of the neighboring 
spin that is pointing in the z-direction) . This then gener- 
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ates spin waves that contain both longitudinal {{S^ {t))} 
and transverse {{{S^ (*))}> {{^n (0)}) components. 

The space-time diagram of {S^ (t)) looks very similar 
as Fig. 2 and therefore we do not show it. Now, we in- 
vestigate the propagation of the transverse spin waves 
by considering one of the two components (the actual 
choice is irrelevant). In Fig. 3, we present results for the 
time evolution of {S^ (t)) for n = 1 to A''. In the Ising- 
like Heisenberg chain without a DW, the transverse spin 
waves propagate in the same manner as the longitudinal 
waves, (compare Fig. 2(a) and Fig. 3(a)). However, from 
Fig. 3(b,c) it is clear that the transverse waves do not 
propagate through the DW structure but are reflected 
instead. 

For a more quantitative study of the interaction of 
DW(s) and spin waves in quantum spin chains, we ana- 
lyze in detail, the time evolution of the right-most spin. 
In Fig. 4 we plot {Sff {t)) as a function of time for the 
six cases depicted in Figs. 2 (Fig. 4(a)) and 3 (Fig. 4(b)). 

From Fig. 4, we conclude that the propagation of lon- 
gitudinal spin waves in the two cases is essentially the 
same, except for the amplitude. Rotating the first spin 
by 7r/2 (instead of tt in the case of the spin flip) about 
the X or y axis generates waves of which the amplitude 
of the longitudinal component at the site N is half of 
that of the spin-flip case. Using Eq. (3) and the fact that 
J25{x) has a first maximum at a; w 27.4, we find that 
(S'2g(t)) has a minimum at t/r « 144. This value is in 
agreement with the time at which the numerical solution 
for the N = 26 chain exhibits a first dip (see black (solid) 
line in Fig. 4(a)). A first conclusion from this analysis is 
that the qualitative aspects of the interaction of the lon- 
gitudinal spin wave excitation and the DW(s) does not 
depend on the transverse components of the spin wave. 

Fig. 4 also clearly shows that the presence of a DW in- 
creases the speed at which the excitation travels through 
the DW. Comparing the curves for the system without 
DW, one DW, and two DWs, we conclude that the solid 
curve lags behind with respect to the dashed curve, and 
the dashed curve lags behind with respect to the dot- 
ted curve. Thus, the longitudinal component of the spin 
wave excitation is shifted forward as it passes a DW. 

Fig. 5 shows the time evolution of (5^ {t)). In contrast 
to the longitudinal component (see Fig. 4), the maximum 
amplitude of the transverse signals strongly depend on 
the presence of DW(s) in the system (note the differ- 
ence in scale between Fig. 5(a) and Fig. 5(b)). Thus, in 
the quantum system, the reflection of the transverse spin 
wave excitation is significantly larger than the reflection 
of the longitudinal component. 

4. Discussion and Summary 

Finally, we point out the difference between the con- 
tinuous model for mesoscale magnetic systems and the 
present lattice model. In the former, DWs exist as rota- 
tion of the spins according to a soliton structure, while 
in the microscopic quantum system, there is no struc- 



ture in the transverse spin component and a DW is de- 
fined as a dynamically stable structure of the longitu- 
dinal components. We found that such DWs exists for 
A > 1 whereas for A < 1 they are unstable. We also 
studied spin wave propagation and found that the longi- 
tudinal components of the spin wave speed up when they 
cross a DW. The transverse components of the spin wave 
are almost totally reflected by the DW, but this charac- 
teristic feature of the microscopic quantum chain is not 
found in mesoscopic magnetic system, where the trans- 
verse components crosse a DW without reflection and 
with a phase shift of 7r/2.^ It should be noted that the 
system described by the LLG equation is fundamentally 
difi^erent from the system that we consider in this Let- 
ter. The former treats the magnetic system in the meso- 
scopic regime as a classical, continuous medium, whereas 
the present study treats the magnetic system as a micro- 
scopic, quantum mechanical system. Which of these two 
approaches is the most suitable description obviously de- 
pends on the specific material. The change of behavior 
from mesoscopic to microscopic may become important 
as bottom-up chemical synthesis is providing new ways 
for further down-sizing of the magnets. 
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